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Abstract 

Wigner and Husimi transforms have long been used for the phase- 
space reformulation of Schrodinger-type equations, and the study of the 
corresponding semiclassical limits. Most of the existing results provide 
approximations in appropriate weak topologies. In this work we are con- 
cerned with semiclassical limits in the strong topology, i.e. approxima- 
tion of Wigner functions by solutions of the Liouville equation in I? and 
Sobolev norms. The results obtained improve the state of the art, and 
highlight the role of potential regularity, especially through the regularity 
of the Wigner equation. It must be mentioned that the strong conver- 
gence can be shown up to 0(log-) time-scales, which is well known to be, 
in general, the limit of validity of semiclassical asymptotics. 
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1 Introduction and main results 

Consider a wavefvmction u £ (x,t) satisfying the Schrodinger equation, 

u £ (t = 0) = u%{x). 

The Wigner transform (WT) W e (x, k, t) of the wavefunction u is defined in 
the standard semiclassical scaling, 

W{x,k,t)= J e- 2myk u £ {x + e^,t)u £ {x-e^,t)dy. (2) 

In that case W corresponds to a pure state. A Wigner function W need not be 
the Wigner transform of a wavefunction u; when working with mixed states the 
Wigner function can be 

W £ (x,k,t)= J e - 2 ™y k p £ (x + e^x-e^t)dy (3) 

i/eR" 

for the density matrix p £ {x,y) 1 which satisfies the von Neumann- Heisenberg 
equation 

ie P £ t + 4(A K - A y )p £ - (V(x) - V{y))p £ = 0, 

(4) 

p £ (t = 0)=p £ o . 

(pQ corresponds to a pure state if there is a Uq(x) £ L 2 such that Pq(x, y) — 
ul{x)u £ {y)). 

Now denote by $ the e-dependent smoothing operator 

*:/(*, fc) (^) " / e~~ L^^ + ^T J /(*', ( 5 ) 

a x and are e-independent parameters that can be used for fine tuning. As a 
matter of notation we will use / = $/. 

Whether we deal with pure or mixed states, the smoothed Wigner transform 
(SWT) is defined as 

W = $W. (6) 
The WT and SWT satisfy equations associated with ([T]), As was shown 
in [3], the equation for the SWT W e (x, k,t) is 



d t W £ + (2irk ■ 3 X + e -^0 x • d fe ) W £ - 



:Re 



i J e 2 ™ Sx ~"r< s V(S)W £ (x + H2k* , fc _ f )dS 
W(t = 0) = Wfi. 



(7) 
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In accordance with the meaning of a x , a k , the well-known equation for W e can 
be recovered by setting a x = Cfc = in ([7]). 

It is well known that given a family of problems as in (fl}, studying the 
limit e -> presents difficulties, both in terms of analysis (e.g. the appearence 
of caustics; also worth mentioning that typically u e doesn't have a meaningful 
limit in e) as well as computation (the complexity of solving numerically problem 
(JTJ) blows up as e -» 0). That's the reason why several asymptotic techniques 
have been developed to treat semiclassical (also called high-frequency, or short- 
wavelength in some contexts) problems - one approach being based on Wigner 
measures (WMs). 

WMs are a very powerful and elegant tool, which under minimal assumptions 
in many cases allows a simple and elegant description of the semiclassical limit. 
On a mathematical level, the idea is that (in an appropriate sense) limW £ is 

a natural object (even when limu e is not), and evolves in time under simple 

asymptotic dynamics. 

A very simple way to state the core result is that, if the potential V is smooth 
enough then (up to the extraction of a subsequence) 

W e {t) W°(t) 

in an appropriate weak-* sense, where 

d t W° + 2nk ■ d k W° - ±d x V ■ d k W° = 0, 

W°(t = 0) = MmWjj, 

see e.g. [12]. The Liouville equation can then be solved with the method of 
characteristics, i.e. is effectively reduced to an ODE problem. 

WMs have been also used to derive asymptotic models in problems featuring 
systems of equations, nonlincaritics, stochastic or periodic coefficients, inverse 
problems etc; see for example [121 [3 [HI HOI H3J HH HH EE [ID] and the 
references therein. However, there are intrinsic limitations to the WM approach, 
that have only recently started to attract systematic attention - see e.g. [B], 
section 11.2.31 

The motivation of the present work is threefold, and can be summarized as 
follows: 

On a technical-mathematical level, we show results that improve upon the 
state of the art on strong-topology semiclassical limits. See section 11.2.11 for 
more details. It must be noted that the application (and extension) of the 
machinery developed in [3] is crucial to the proofs here. Moreover, all the error 
estimates are constructive, most of the constants are computed explicitly - in 
fact all of them are computable, including in particular the time dependence. 

On a numerical level, our results justify the use of the SWT for fast coarse- 
scale simulations; see e.g. [TJH], and section [1.2.21 
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Finally, on a more qualitative level, the ideas developed here can be used as 
an intermediate step to work out problems in which the semiclassical limit is 
not effectively known; see section [I.2.3l 



1.1 Statement of the assumptions 

Let r € N U {0} be fixed. The main assumption will be 
Assumption Al(r) There is a Mo > such that 

J \V(S)\(1 + \S\ r + 2 )dS = M . (8) 

Now denote the auxiliary potentials V\(x), V\{x) to be 

^( x ) = (f) f Je-^ lx - x ' l2 V(x>)dx>, 

M*)=(4) f je-^ x - X ' ? V,{x')dx' 
The parameter rj is set to be 



In addition, Vrn = 0, 1, .., r 



(9) 



V = \ 2" mCLX V U ~ ^ ><Tx ' ~~ + 1 ) £ 

(see equation I1101\) and above for the rationale behind this scaling). 
There are C,D,R> 0, 9 G (r + 1, r + 2) such that 

VmeN,m<0+l, : |A| < m ||fl£^i(aO||L-> «S C, (10) 
Vm G N,m > (9 + 1, : |A| < m => ||d^Vi(a;)||L°° < 



(11) 



(l-e-<^' s )\S\ m+1 V(S)\\ LHMn) ^C(e^+e) (12) 



There is a somewhat stronger version of this assumption which is much 
simpler: 

Assumption Al'(r) There is a Mq > such that 

J\V(S)\dS = M , (13) 
and 3R, D > 0, 9 E (r + 1, r + 2) suc/i </ia< 

V|5|>i? : < (14) 



We mentioned that Al'(r) Al(r). Indeed, 
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Lemma 1.1 Equations f73J) and (7J]) imp/?/ 0), Q7J), fUJ, UM>- 
Proof: The proof is given at the beginning of section 14.41 

In some cases we will use 
Assumption A2 Al(r=0) holds, and moreover ||Wq||#i = o(e~5 \\Wq \\l 2 )- 

Remarks: 

• The reason we use the more complicated (but weaker) conditions (flO|) . (fTTj) . 
([12")) instead of (fl"4"|) is that there are relevant examples not covered by it; a simple 
example is given by a lower-dimensional potential: consider V*{x\, x n ) — 
V* (x%, Xd) with d < n; if 14 satisfies 

-(d+x+fl) 

/ f ' \ 

V|S| >R : 

it is easy to check that Al follows - although in general (fl"4"|) does not. 

• Following assumption Al(r=0), V(x) £ L°° . This suffices to imply that the 
Schrodinger operators — A A + V(x), — kA + Vi(x), — \A + V\{x) are essentially 
self-adjoint on L 2 (M. n ). This is something that we will use, e.g. in applying 
theorem 12.11 (and the reason why we introduce Vi here). 

We will not comment on this explicitly any more. 

• It is clear that if Al(r) is true, then Vr' <E [0, r] n Z Al(r') is also true. 
1.2 Formulation of the main results 

Theorem 1.2 (L 2 semiclassical asymptotics) Consider a potential V(x) and 
an initial Wigner function Wq , such that (A2) is satisfied. 

Denote by W E (x, k, t) the solution of the corresponding Wigner equation with 
initial data Wq as in (A2), W e (t) = $>W e (t) the corresponding SWT, p(x,k,t) 
the solution of the Liouville equation 

d t p + 2nk ■ d k p - ±d x V ■ d k p = 0, 

(15) 

p(t = 0) = W$, 
and by p\(x,k,t) the solution of 

d t pl + 2nk ■ d k p\ - ±dj/x ■ d k pl = 0, 

( 16 ) 

fi(t = 0) = W§. 

Then there is an 0(1) explicit constant C, depending only on n, a x ,ak } Mo 
such that, for all t € [0, T], m G N, 

\Mt)-w%t)\\ L * < 

(17) 
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\\pl(t)-W e (t)\\ L2 n rzrj, Tn-K max{l,\\V{S)\S\ 2 \ \ r t \ Wo\\jl}_ _ 

\\w S \\ L2 ^W^e \\w S \\ L .- 
= o (r e T™rnax{i,\\v(s)\s\ 2 \\ Ll }^ , 

\\p(t) - W s {t)\\ L * ^ 
s: C(Ve(T + + Te\\W^\\ L 2)e Tn ^ max ^^^ s \ 2 W^\ 

\\p(t)-W e (t)\\ L2 r r-,-. T x Tmrmax{l,\\V(S)\S\ 2 \\ rl } IH^ikl - 

\\w S \\ L 2 ^ oveU + J je \\wg\\ L 2 - 

(20) 

T yTn^max{h\\V(S)\S\ 2 \\ Ll }' 



(18) 



(19) 



(<> 



Remarks: 



• Observe that the result can be extended beyond time-scales T = 0(1), up to 
T = 0(—log(e)). As is well known, semiclassical asymptotics in general cannot 
be extended for longer time-scales [3 [17] . It is worth noting that we can go to 
such time scales with Al(r=0), which amounts to 2 derivatives on the potential, 
as opposed to analytic potentials (in [4]) or C°° potentials (in [5]). 

• Similar approximations hold e.g. between W e (t) « pf (£) , etc; see the proof. 

• One can see in the proof that all constants are (or can easily be) computed 
explicitly. Although of course they will be pessimistic, given a problem (i.e. 
potential and initial data) they provide one with an a priori estimation of all 
the errors involved. This is true for the following theorems as well. 

• Pure states, i.e. Wq = W £ [u £ ], typically give ir^ii^ = 0(e~t ) - see lemma 
IB. II Therefore this result does not yield a significant estimate for pure states 
(i.e. the relative error bound obtained is not o(l)). See corollary 11.61 and the 
remark thereafter for pure states. 

Theorem 1.3 (Negative-index Sobolev spaces) Consider some r £ N. As- 
sume that Al(r) is satisfied, and Wq £ H 1 ^ r . 
Then, for all t £ [0, T] 

\\p £ (t)-W £ mH-^CT\\W^\\ H ^-re TD ^ n \^-e + ei)+ 

+ C{e TD(r,n) + gTpM+DM.n)]^ | | + (21) 

+ (e + ei)CTe T ^ r+1 ^+ D ^\\Wi\\ H -i-r. 

As long as there is a constant C\, independent of e, so that \\W e (t)\\jj-r ^ 
C\, equation i21\) yields a relative error estimate in a straightforward manner. 
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Remark: As long as the assumtpions of [T2] are satisfied, it can be used to 
bound from below ||W E (i)||#-r. Take (j> € S, \\4>\\H r ^ 1; [H] assures that 
(W £ , 0) -> (W°, 4>). So a sufficient condition is that there is a (f>o{t) so that, for 

te[0,T], KW°(*).0o(*))| >2Ci>0. 

If that's true, then it follows automatically that 

3e :Vee (0,e ) : |<W°(t) - W £ (*), O (*)>| < d- 

Of course this is mostly a technical point; we are not aware of any meaningful 
example where \\Wq \\h->- = 0(1), and later ||VF e (£)||#-r = o(l). 

Theorem 1.4 (Positive-index Sobolev spaces) Consider some reN. is- 
sme i/iat Al(r) is satisfied, and Wq £ H +r . 
Then, for all t €E [0, T] 

||^(t) - W*{t)\\ H r <C\\W§\\^ [Te TD ^ n \^-e + ei) + 

(22) 

+Te T[D(r+l,n)+D(r,n)]( £ + £ V 1 ) + e TD(r,n)^ 

Moreover, if in addition 

Wl*- -«.-»), (23) 



||^(t)|| ffr + l 

equation \22]) yields directly a relative error estimate as well. 

Remark: A sufficient condition for equation (|2"5|) is given by ^[j^if^ 1 
o(e~i). Indeed, to see that, apply 

\\W%t)\\ Hr+1 > \\W e (t)\\ L3 => — — L < 



to estimate Mwe ,.| M in the lhs of equation (1231) . 



It is interesting to look at some corollaries: 

Corollary 1.5 (L 2 limit) Assume that (A2) holds, and 3W§ <E L 2 (R 2n ) such 
that 

lim\\W° -W£\\ L 2 =0. 
Denote p°(x,k,t) the solution of 

d tP ° + 2nk ■ d kP ° - ±d x V ■ 8 k p Q = 0, 

(24) 

p°(t = 0) = VK °. 



Then, Vi e [0,T], 



lim||p°(t) - W £ (t)\\ L 2 = 0. 

£->0 
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Corollary 1.6 (Strong limit with concentration) Assume thatBs £ N, W§ £ 
H 1 ^ 3 such that 

lim\\W§ - W E || H - 8 = 0. 

Assume moreover that assumption Al(s) is satisfied, and recall the definition of 
p°(x,k,t), in equation \21$ . 
Then, Vt G [Q,T], 

Hm||^ (i)-W £ (i)|| ff -. =0. 

Remark: Pure states. It is important to note that this result can also be 
applied to pure states: for example if the initial data is a coherent state, then 
(provided the potential is smooth enough) the assumptions of corollary 11.61 are 
satisfied for s = + 2; see lemma HOI 

However, strong aproximation doesn't have to mean strong limit: 

Corollary 1.7 (Strong approximation with concentration) Assume (A2) 
holds, and Wq — 1 Wq (fc L 2 in weak-* sense (as in J12j). 
Then 

Remark: While neither p £ (t), nor W £ (t) have a limit in L 2 as e — > 0, they are 
asymptotically close to each other in L 2 Ve > 0. The point here is that there are 
situations where this might apply while corollary 1 1 . 61 not - for example when the 
convergence to the Wigner measure isn't known in norm, or when the potential 
isn't smooth enough. 

Observe finally that an L 2 approximation is much stronger than and H~ s 
limit. For example, the L 2 approximation implies any oscillations that might 
develop during the quantum evolution (i.e. absent from the initial data) are 
necessarily small in L 2 sense - while of course the iJ~^"f 1+2) ]j m it tells us only 
that eventually any oscillations will cancel out in weak sense. 



Finally there are a couple of somewhat technical points worth mentioning: 

• As long as a x ak ^ 1, the SWT "preserves positivity", i.e. if W e (x,k) is the 
Weyl symbol of a positive operator, then W £ (x, k) ^ 0. (This is well-known, and 
is a property that the WT itself does not have; see lemma lA~2l) . Our asymptotic 
analysis here does not assume anything on the sign of W E (x, k). Still, it should 
be noted that a strong approximation which preserves sign might be of particular 
interest in more difficult (e.g. nonlinear) problems. 

• It is also important to note that the basic limitations come from the break- 
down of Sobolev regularity in pf, W e . Therefore, to understand where there 
is an actual breakdown of the approximation, as opposed to a mere technical 
limitation of our tools, it will be important to get sharper assumptions on the 
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validitiy - and failure of theorems such as 12.21 The weaker assumptions e.g. 
of theorem II .91 is one possible way forward. 



1.2.1 State of the art 

One reason for working out strong topology semiclassical approximations for the 
Wigner function, is that there is an extremely limited literature on the subject. 
Indeed, most results focus on weak-* limits. 

From an analytical point of view, a norm could quantify a rate of convergence 

- in contrast to the weak-* convergence, which depends on the choice of test 
function. This could also be extremely important with respect to computational 
applications (see next section), where it is important to quantify a priori (even 
if on the pessimistic side) the errors involved. From a more qualitative point 
of view, it could be the way to add e-dependent corrections to the limit, thus 
controlling in a more flexible way the information loss. 

As long as we are working with norms, L 2 norms are especially well suited 

- in particular in view of theorem 12.11 The H m generalizations also, in fact, 
make use of theorem 12.11 (L 1 approximations also make sense, in view of the 
classical interpretation. We don't work on that here, but it may be a direction 
worth pursuing in the future). 

To the best of our knowledge [TH] comprises the state of the art for strong 
semiclassical approximations. Other existing works (e.g. [3]) often use regularity 
assumptions stronger than [19] , and, to our knowledge, do not treat H m norms. 
It is therefore natural that our results here should be compared to those of [19] . 

In our notation, Theorem 5.1 of [19] states (or rather implies, since it also 
treats higher order approximations): 
Statement If 



However, there is a misprint in |19j ; the statement of the theorem should 
read 

Theorem 1.8 7/ 




(25) 



then 



\\p(t)-W*(t)\\ L , ^0(e 2 \\W^\\ H ,e 



ct(i+f\v(S)\\s\ 2 dS) 



) • (26) 




(27) 



then 



\\ P {t)-W^t)\\L^o[e 2 




.Ct(l+f\V(S)\\S\ 3 dS) 



) ■ (28) 
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(The problem in Theorem 5.1 originally comes from from a typo in Theorem 
4.1. See section[2] theorem l2.2l for a detailed treatment of the point in question). 

In that light, the use of the SWT in Theorem 11.21 allows one to gain one 
derivative in each of the initial data and the potential. (For example V(S) = 
5— r satisfies (A2), while f \V(S)\ \S\ 3 dS = oo). 

The price one has to pay for these weaker regularity assumptions, is a larger 
size of the error: for example, if > 2 and ||Wq||#2 = O(l), theorem 11.21 gives 
\\p(t) - W 6 (t)\\ L 2 < O(Ve), as opposed to an 0(e 2 ) error from Theorem 5.1 of 

It must also be noted that neither result applies to pure states. [TH] contains 
a result in H~ m developed with pure states in mind, Theorem 5.2; our approach 
here also yields an H~ m approximation result applicable to pure states, corollary 
11.61 which follows from theorem 11.31 (Note that theorem 11.31 also improves on 
the regularity assumptions of Theorem 5.2 in [19]; we don't do the exhaustive 
comparison here). 

1.2.2 Numerical applications 

Part of the motivation for this work has been its potential for numerical ap- 
plications. Indeed, a numerical scheme using the SWT for the fast recovery of 
coarse-scale observables was formulated in [TJ [3] , and proof-of-concept numerical 
experiments showed that it is a competitive option in the semiclassical regime. 
This work serves as a rigorous justification for the method. 

In this regard it is worth highlighting one point in particular: it is well known 
that keeping track of the WM allows one to approximate quadratic observables 
of interest, with famous example 

\u e (x,t)\ 2 « J W°(x,k,t)dk. 

In addition, WMs don't break down - or even require any special treatment - 
when caustic appear. However, although a WM stays well defined on phase- 
space, its marginals need not exist, especially when a caustic appears. That is, 
while 

Ve>0 \u £ (x,t)\ 2 = J W £ [u £ (t)}(x,k,t)dk, 

J W°{x, k, t)dk may not be well defined, even if lim \u £ (t)\ 2 ifl See [9] for more 

details. This means that, if we are inetersted in an approximation of \u 6 {x, t)\ 2 , 
and not the WM W° itself, the WM technique does, in fact, have problems on 
caustics. These problems can be overcome using the SWT. 

This raises a more general question: can we build e-dependent corrections 
of some sort into a phase-space asymptotic technique - in a computable way? 

1 lim \u e (t)\ 2 = +00 is still considered well defined; J W°(x, k, t)dk may give rise to expres- 
sions that have no value, finite or otherwise 
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Recall that W° does not depend on e, i.e. any WM-based technique cannot 
offer any kind of e -dependent information. 

This is the motivation for constructing, in theorem 11.21 the approxima- 
tion p\ « W e , in addition to p « W e . From a mathematical point of view, 
p(t) = W £ (4>t(x,k),t = 0) is a good approximate solution. From a numerical 
point of view however, the oscillations that W £ (x, k, t = 0) in general contains 
could make even its propagation along characteristics as difficult as the direct 
solution of the original problem ([!} - keep in mind that this is in twice the space 
dimensions than (J]). On the other hand, the oscillations of W§ are controlled 
by the smoothing, making the propagation p%(t) = W e (<fit(x, k),t — 0) signifi- 
cantly easier. (The fact the the flow comes from a smoothed potential also can 
simplify a number of things, e.g. in setting up a solver for the trajectories). So 
finally smoothing gives a practical way to incorporate e-dependent information 
in a phase-space asymptotic technique. 

One should keep in mind that in general the complexity of simulating semi- 
classical problems explodes as e -> 0, so apart from documenting the convergence 
and limitations of any mumerical scheme, a key question is to understand the 
rate at which complexity (computation time, memory used) grow in s. See [UH] 
for more details on the use of the SWT in semiclassical simuation, including in 
particular a study of computational behavior in 1-dimensional problems, of the 
approximation of observables (with main example \u e (x, i)| 2 ), and examples of 
calibrating the smoothing with respect to intrinsic scales of the initial data. 

Finally, the fact that all the errors in this work are computed explicitly (in 
contrast e.g. with weak-* results, where the errors depend, in a non-explicit way, 
on the test function) can be very helpful in setting up numerical simulations. 

1.2.3 Further work 

To conclude this long introduction, some words are in order for the more indirect 
motivations and implications. It has become more and more clear that there are 
problems where WMs, in their more straightfoward application, fundamentally 
fail. (For example, when there can be constructed two (e-families of) problems 
whose WMs are identical at one point in time, but not at a later one). 

Examples include nonlinear problems, systems with crossing eigenvalues etc 
- see [5J for a survey and a plethora of references. Another example, triggered 
by low smoothness in the potential, is presented already in [T^], Rem. IV 3. 
Loosely speaking, in all these cases the information lost in the limit to W° is 
"too much" ; it is worth wondering if some modified, augmented version can still 
be made to work. In some cases this has already been done; see e.g. [BJ. 

The following is a version of theorem 11.21 under weaker assumptions (one can 
easily check that the proof still applies) : 

Theorem 1.9 Assume that (Al) holds for some 6 £ (0, +oo)\N (i.e. 9 possibly 
smaller than 1 ), and in addition that there is a constant M\ such that for all 



11 



te[o,T] 



n 



\ = o((e-S+e-i)\\W*(t)\\ L ,) 



(29) 



See equation J77| ) and thereafter for the precise definition ofW((t). 
Then there is an 0(1) explicit constant Cq such that for t £ [0,T] 



P \{t) - W e (t)\\ L 2 < T Mi Co (e* + Ve) , 



(30) 



and 



M(t)~W%t)\\ L 2 

\\W°{t)\\L* 



0(1). 



(31) 



By modifying this kind of result (with nontrivial differences), we are able 
to prove a result applicable to problems with potentials satisfying (Al), (A2) 
for some appropriate 9 6 (0, 1). It must be noted that 9 < 1 vs. 9 > 1 is a 
nontrivial dichotomy: 9 £ (0, 1) includes C 1 \C 1 ' 1 potentials, for which there are 
some results, see [12], but they show convergence to an ill-posed semiclassical 
limit. In other words, there are situations with < 1 where we can show 
that (an appropriate version of) condition (|2T))) holds, which finally allows the 
regularization of the known, but ill-posed, limit problem [work in preparation]. 
A lot of the ideas and tools that had to be developed for that are in fact in this 
paper. 

That is, apart from an end in themselves, the results shown here (including 
in particular the regularity results, section [5]) are necessary in studying other, 
more challenging problems as well. 



The rest of the paper is organised as follows; below we define some notation 
we will use; in section [2] we will present a family of results concerning the 
regularity of quantum (Wigner) and classical (Liouville) equations, trying to 
be as thorough as possible. The proofs can be found in section [3] Section 
0] is devoted to the proofs of the theorems 11.21 11.31 and 11.41 The corollaries 
stated earlier are straightforward applcations of these theorems, and theorem 
ll.9l follows readily by retracing the proof of theorem II .21 

1.3 Notations 

The Fourier transform is defined as 




f(x)dx. 



(32) 



12 



Inversion is given by 

/(*)= ^xi* [/(*)] = / e 2 « ik *f(x)dx, (33) 

^ [^a\ [/(«)]] = ^ ^ [/(«)]] = /(*)■ (34) 

The Sobolev norm of order m on phase-space will be defined as follows: 

II/IIh~(R>») = E \\9 a Af\\ m ^), (35) 

|o|+|6|<m 

where of course a and b are multi-indices of length n each. More generally we 
have 

||/||*»(r«) = E H^/IU 2 ^)- (36) 

|a|^m 

H~ m will be the dual of H m . This definition is equivalent to the more usual 

ones, e.g. based on ||/(ar)(l + M) m || L 2, or / E ii^Tiliw^r 

Vlo|<"» 

A more general family of Sobolev- type spaces that we will use are X m ' p , 
with norm defined as 

\\f(x,k)\\ xm ,» = \\(\X\ 2 + \K\ 2 )%f(X,K)\\ L *. (37) 
Moreover, denote for future reference 



Finally, we must mention an abuse of notation that we will do: sometimes 
we will suppress the e dependence of certain functions for economy, using e.g. 
W(t) instead of W e (t), px{t) instead of p\{t) etc. 

2 Quantum and classical regularity results 

The following theorem is a direct consequence of Theorem 2.1 of [13] : 

Theorem 2.1 (L 2 regularity of the Wigner equation) If the Schrddinger 
operator — + V(x) is essentially self-adjoint on L 2 (W l ), then the correspond- 
ing Wigner equation preserves the L 2 norm, i.e. for Wq € L 2 (M. 2n ) there is a 
unique solution of 



d t W e + 2irk ■ d x W e + ^Re 



i J e 27riSx V(S)W e (x, k - ^)dS 



W e {t = 0) = w e , 

and HW^tJHjapp.) = \\W£\\ L 2 (M 2 n) V* € R. 



°' (38) 



13 



The following estimation is quoted from [19J . where it appears as Theorem 
4.1 (with a different but equivalent definition for H m ): 

Statement (H m regularity of the Wigner equation) For each m 6 N 
denote 

C m = J\V(S)\ \S\ m dX. 

R n 

Denote by W E (t) the solution of the Wigner equation H38\) , with initial data 
W e (t = 0) = Wq . Then, if C m < oo, there is a constant C independent of e 
such that 

\\W s (t)\\ Hm <e ct |ra ff ™. (39) 



However it contains a typo. We prove a weaker result: 

Theorem 2.2 (H m regularity of the Wigner equation) Denote by W £ (t) 
the solution of the Wigner equation i38\) , with initial data W e (t — 0) = Wq . 
Then, for each m G N, the following estimate holds for time-scales t £ [0,T], 
T = 0(l):if 

r \V(S)\ (|5| + |5| m+1 )^<oo, 



then there is a constant D — D(m, n) such that 

\\W £ (t)\\ Hm ^e tD ^\\W^\\ Hm . (40) 
Remark: See equation for an estimate of the constant D(m,n). 



The same strategy as in the proof of theorem 12.21 works for the Liouville 
equation as well: 

Theorem 2.3 (H m regularity for the Liouville equation) Denote by p(t) 
the solution of the Liouville equation 

p t + 2TTk-d x p-^-d x V(x)-d kP = (41) 

with initial data p(t = 0) = po- 
Assume that 

wv(s)\sr +i \\ L .<^ 

and that, Vt € [0, T] 

| (d^p{t),k ■ d x d^p{t)) | < +oo, 

(42) 

| (d£d* P (t), d x v ■ d k d£d* P {t)) | < +oo. 
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Then the following estimate holds for t € [0,T]: there is a constant C > 
such that 

\\p(t)\\ Hm ^e tc \\ Po \\ Hm (43) 

A sufficient condition for equation \J^2^ to hold is the following: 
The flow associated with J^_?[ ) will be assumed to be complete, i.e. any bounded 
ball is mapped within some bounded ball for all times ( trajectories do not reach 
infinity in finite time). Moreover that po is of compact support and, for every 
multi-index a with \a\ ^ m + 2 

d a x V(x) € L&CR")- (44) 

Another sufficient condition is that for every multi-index a with \a\ ^ m + 2 

d a x V{x) £ L°°(R n ). (45) 

Remark: Our sufficient condition invoives one more derivative in the potential 
(up to order m + 2), but the estimate itself involves only derivatives of order 
TO + 1. This will be important for example in situations where the derivatives 
exist but become asymptotically large (e.g. if V was a mollified version of a 
non-smooth function). 

There are some more regularity results for the Wigner equation in section 
13.31 We don't discuss theme explicitly because we don't use them directly in 
the proof of our main results here. However it must be noted that in other 
situations they might be very useful. 

3 Proof of regularity results 

3 . 1 Proof of Theorem B 

The Fourier transform of the Wigner equation (|38p is 

W t - 2ttX ■ d K W + 2 J V(S)W(X - S, K ) sm ^ *0 dS = ( 46 ) 

For future reference denote by U(t) the propagator of this equation. Theorem 
12. II essentially tells us that, under the self-adjointness assumption, 

VtGR \\U(t)\\ L 2 {M 2 n ^ L 2 (R 2 n) =l. (47) 

Denote 

v A , B (X, K) = X A K B W(X, K) (48) 
in the usual multi-index notation. By elementary computations, it follows that 

X A K B X ■ d K W = X ■ d K (v A}B ) - £ B : r x .. j. , (m 

Bj>0 V ' 
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and 

X A K B J V(S)W(X - S, K ) sm ^f- K) dS = 
= J V{S)v a .b{X - S, K) sln ^f K) dS+ 

+ E ft (t) JV(S)S l v A -iMX - S,K) ainln f K) dS. 

0<l^Ai=l * 

To see that, it suffices to observe that 

A 



(50) 



x A = {s + {x - s)) A = E (1)s l (x - s) 



A-l 



1=0 



=> x A = (x-s) A + E H)s l (x - s) A - 1 

10L(A) 



where for brevity we use the notation 



and 

h(A) -.= {i e n n \i s? A, \l\ > 0}}. (51) 

n 

(Remark on notation: I ^ A Vi € {l,...,n} : ^ Ai, and |Z| = E k)- We 

will use these notations freely in the sequel. 

Two things should be noted: obviously if Ai = then there is no contribution 
from the i coordinate; and / € h(A) -> < \l\ ^ \A\. 

So multiplying equation (|46|) with X^if 5 yields 

|va,b ~ ■ d K v A . B + 2 / V{S)v A , B (X - S, K y^f K ) dS = 



Bj>o (52) 

-2 E OfV(S)S l v A ^ B (X - S,K) sm{ *f K) dS, 
leh(A) 

which, for any m £ N, is a closed system for v At B with \A + B\ ^ m. (For 
A = 0, we just forget the last term). 
Moreover, observe that 

j t \\vAMt)\\h = 2\\vA, B (t)\\v J t \\VAAt)\W (53) 



1G 



while at the same time, using equation 



jt\\vA,B(t)\\ 2 L 2 = (v A ,B,d t v A<B ) + (dtv A ,B,v A , B ) -- 

= 2Re [(v AiB (t),2irX ■^ K v A , B )- 
-(v A Mt),2 I V(S)va,b(X - S, K) sm ^f K) dS)- 

+ (v A . B ,2lT BjVA+e^B-ej)- 

Bj>a 



(va,b, E tf)2fV(S)S l v A -iM x - s > K ) 



sin(-KeS-K) 



dS) 



Assuming for a moment that the integrals exist, we observe that 
Re [(v A , B (t), 2irX ■ d K v A , B )] = 0, 



and 



Re 



(v A , B (t),2 I V(S)v AB (X-S,K) Sm(7T£S ' K) dS) 



= 



by anti-symmetry. 

Equation ([53]) is completely obvious; for ([55)) observe that 

Re [j f(X, K)2 J V(S)f(X - S, K) sm{n f K) dSdXdK 

= JJ(X,K )V(S)f(X - S,K) sm{7T f- K) dSdXdK + 
+ J }{X, K)V(S)f(X-S,K) 



sin(7rsS- K) 



dSdXdK = 



= ff(X,K)V(S)f_(X-S,K) 
+ Jf(X,K)V(-S)f(X-S,K) 



e 

sin(7T£S-K) 



dSdXdKi 
dSdXdK 



(54) 



(55) 
(56) 



= j f(X,K)V(S)f(X - S,K) sm(v f- K) dSdXdK- 
- J f(X, K)V{S)J(X + S, K) sm ^f- K) dSdXdK = 

= fJ(X,K)V{S)f(X - S,K) sm(7T f- K) dSdXdK- 
- J f(X - S, K)V(S)J{X, K) sm{n£ J- K) dSdXdK = 

Of course for these considerations to be valid the integrals must exist; they 
do, thanks to the following lemma: 

Lemma 3.1 (H" 1 regularity of the Wigner equation) Denote byW(t) the 
solution of the Wigner equation like earlier. Assume that 



\V(S)\ (\S\ + \S\ m )dS < oo, \m\\ H m < e 



(57) 
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for some c > 0. 

Then, for each m 6 N, T > 0, i/iere exisi constants C, /i > suc/i i/iat 

||W e (t)||j*« < Ce"", (58) 

and 

V\A + B\ <m-l, i€{l,...,n}, \\hd£d% W(t)\\ L 2 < Ce - **. (59) 

(For our purposes it is obvious the explicit computation of C, /it is not im- 
portant - it is possible, as one can readily see in the proof of lemma 13. 1[ in 
section |3~51 Observe moreover that the requirement V(5)|S'| 6 L 1 , which allows 
for weaker possible blowup at than V A (S , )|S , |™ 1+1 G L 1 , originally comes from 
lemma 13. 1|) . 

Now, by combining equations (|53p and ([51)1 and discarding the terms ([55)) . 
that vanish identically, it follows tha10 

Tt\\ V A,B(-t)\\L* E ^ll^+e„S-e,(t)||^ + 

S 3 >0 

(60) 

n 

E (f)Elins)s i+e *IMK-*,iW*)IU- 

Now we sum up all the equations for |A| + \B\ ^ m to get 

i E IK-bWII^ < 

\A+B\^m 

(61) 

<n m+1 7rm!-maa ; {l ) ||y(5)|5r+ 1 || il } E IKsWIU- 

\A+B\^m 

With the help of Gronwall's lemma the result follows. For brevity, in the 
sequel we will denote 

D(m,n) := (1 + n m+1 m\)ir ■ max{l,\\V(S){\S\ + |^r +1 )|| L i}- (62) 

Remark: For m = 1 it is very easy to see that we can have a slightly better 
constant, namely: 

Tt E \\va,b®\\v < 

|A+S|<1 

(63) 

^rmmax{l,\\V(S)\S\ 2 \\ L i} E IKsWII^, 
i.e. D(1,ti) = nnmax{l, \\V(S)\S\ 2 \\ L i}. 

2 We also used the obvious estimate \v(S)S l v A _ liB (X - S, K) ain ^eS-K) | ^ 
\v[S)S l v A ^ l)B (X - S,K)nS ■ k\ in equation (S3). 
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3.2 Proof of Theorem [2731 

Proof: The proof is essentially the same as for theorem 12.21 we have 



instead of 



7T J V(S)p(X-S, K)S-KdS 
Jv(S)W(X-S,K) Sm{7r£ £ S - K) dS 



in equation (|46p. and then the proof follows along similar lines. The counterparts 
of the terms equations (|55l) , ([55| arc 

Re [{d£d*p,k>d x d£d%p)], 

(64) 

Re [Wft^^ffi]. 

It is clear that, if the integrals exist, these terms vanish identically by anti- 
symmetry. 

Now we will work out our sufficient condition for the existence of the integrals 
of equation (|64|) : it follows easily, using the method of characteristics, that if 

sup |9£V(x)| < oo, 
\a\ ^ m + 2 
x e R" 

and the flow is complete (see the statement of the theorem), then, for initial 
data of compact support such that 

d s d£d?p(0), d x d£d£p(0), kd x d^p(0) e L 2 , (65) 

we have Vf > 

d x d£d* P (t), d x d^ P (t), kd x d*d*p(t) e l\ (66) 

and hence the integrals in question exist. 

The second sufficient condition follows clearly. 
The proof is complete. 

3.3 Auxiliary lemmata 

It is more convenient to prove lemma IXT1 first for m = 1 and then for m G N: 

Lemma 3.2 (H* regularity of the Wigner equation) Denote byW £ (t) the 
solution of the Wigner equation 138\) , with initial data W e (t = 0) = Wq . Then 
the following estimate holds for t £ [0,T]: if 

\V(S)\ \S\dS < oo, 
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then we have 

\\sd Xi W e {t)\\ L2 «S \\sd Xi W$\\ L 2 + 2t\\V(S)\S\\\ L x\\W§\\ L * 

\\ed ki W e (t)\\ L 2 ^ \\ed ki W^\\ L 2 + 

+2nt (\\d Xi W$\\ L2 +2±\\V{S)\S\ \\lA\W8\\v>) , 

\\ekiW s (t)\\ L 2 < \\ek i W£\\ L2 +2irt\\V(S)\S\\\ L i\\W£\\ L2 . 



Proof of lemma I3.2t 

Denote Ui(X,K) = XiW(X, K), Vi(X, K) = KiW(X,K), z t = d Ki W{X,K). 
Then one checks that equation (|46p implies 

d tUi - 2ttX ■ d K Ui + 2jV(S) Ul {X - S,K) sm(7r f- K) dS = 
= -2 J V(S)SiW(X - S, K) sm ^f K ) dS, 



(68) 



8m - 2ttX ■ d K Vi + 2 J V(S)vi(X - S, K) sini ™ s - K) dS = 

= 2-KUi, 

dm - 2ttX ■ K z l + 2 J V(S)zi(X - S, K ) sm{ *f- K) dS = 
= -2tt J V(S)SiW(X - S, K) c °^s-K) e dS _ 

Making use of the second part of lemma l4~2l (equation (|15ip ), it follows that 
an equivalent reformulation of equation (|68[) in terms of the propagator U(t) is 

Ui(t) = U{t)ui(0) + 2 J U(t -t)J V(S)SiW(X - S, K, r) sin ^f K ^ dSdr, 

T=0 

Vi(t) = U(t)vi(0)- J U(t-T)2Tr Ui (r)dr, 

T=0 

Zi (t) = U(t)zi(Q)+2iT J U(t - t) j V(S)SiW(X - S, K, r) cos ^f- R ) dSdr. 

(69) 

Now, making use of theorem l2.1[ the obvious bounds \sin(ireSK)\, cos(ireSK)\ ^ 
1), and lemma ETT1 we readily deduce from equation (IBTJl) that 



\\ui(t)\\ L 2 «S || Ui (0)|| L 2+2|||^(S)|5||| L i||W £ | 



\\Vi(t)\\ L 2 < 11^(0)1^2+2^ SUp ||tti(T)|| ia < 

re(o,t) (70) 
< ||«i(0)|| i2 +2nt(\\ui(a)\\ L i+2i:\\V(S)\S\\\ L i\\W$\\ L 2 



\\zi(t)\\v < + *?\\V(S)\S\ \\lA\W§\\v>. 
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The proof is complete. 



Now we are ready for the general case. 
Proof of lemma I3.lt 

System ([52]) can cquivalently be recast as 

VA, B (t) = U(t)v A , B (0) - J U(t-T)\2TT E B j V A+ejjB - ej (T) + 

(71) 

+2 E (f ) / V(S)S l v A -i,B (X - S, K, r) sin ^f - K) ds] dr. 

ZGL(A) J 

The point to repeating the idea of the previous proof, is to go through the 
multiindices (A, B) so that we always have on the right-hand-side a quantity 
that we have already estimated. 

This is easy to check; for example let us prove the result for m > 1 assuming 
it holds for m — 1. First work with all the indices such that \A\ + \B\ = m, 
\B\ = 0. It is clear that, for those terms equation ([71]) is just 

?U,o(t) = U(t)v Afi (0)+ 

(72) 

+2 / E (f ) J V(S)S l v A _ li0 (X - S, K, r) sin ^f ' K) dSdr. 

But \l\ > (recall the definition of L(A), equation (J5TJ) ) means that all the 
terms of the form v A -i t o in the right-hand-side are of order at most m — 1, 
and therefore bounded by assumption. By using once again lemma 14.11 and 
integrating in time, we get 

\\v Afi (t)\\ L 2 «c ||«a ) o(*)||l»+ 
+2| E tf)\\V(S)S l \\ L i sup \\v A . lfi (r)\\ L 2 < 

ZeL(A) rE(0,t) 

(73) 

^ ||«a,o(<)||l«+ 

+2|||V(S)|Sr|| Ll E (f) *«P ||«A-i,o(r)|U.. 
ieh(A) re(o,i) 

Now we can allow \A\ + \B\ = m, \B\ — 1, and observe that the rhs contains 
only terms of the form \A\ + \B\ = m, \B\ = 0, which we just estimated. To 
conclude we just proceed inductively until \B\ = m. 

The second inequality in (|58|) follows similarly, by a straightforward adap- 
tation of the previous proof. 

The result now follows. 
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Lemma 3.3 (H m regularity of the Wigner equation) Under the assump- 
tions of theorem \2.2l 



\\W s (t)\\ H - m ^e tD ^\\W^\\ H - m . 
Proof of lemma l3.3t The result follows directly by noting that 



sup \(U(t)f,g)\= sup \(f,U(-t)g)\^\\U(-t)\\ H ^„ m \ lfllH -r. 



Finally the following regularity result can be shown with minimal assump- 
tions on the potential - although apparently it cannot be generalized to higher 
order derivatives: 

Lemma 3.4 Assume there is aC* > such that V{x) > —C*, and the Schrodinger 
operator ^-A — V(x) is self adjoint. 

Moreover, consider a mixed state as follow^: denote A 6 A an appropriate 
index set, and 



P%{x,y) = J ul(x)ui(y)dv(X), (74) 
where \\u\(t ~ 0)\\l2 ~ 1, and 

n 

\(V(x)ul(t = 0),ul(t = 0))| +J2 (\\d Xi nl(t = 0)\\v + = 0)\\ L .) < C x . 

i=i 

Moreover, assume (1 + C\) G -L^d/^A), A). (In particular no assumption is 
made for the behaviour of C\ in e). 

Naturally the relation between p e and the Wigner function is as in equation 

W- 

Then there are C , c > such that for all i = 1, n, t £ [0, T], 

\\d Xi W"{t)\\ L ^C'e- c , 

\\d ki W%t)\\ L2 ^C'(l + T)e- C . 

Proof: First of all, if u\(x,t) is the solution of ([T]) with initial data u\(x), it is 
straightforward to check that 



p E (x,y,t)=J u £ x {x,t)ul{y,t)dn{\), 
i.e. the measure dfi(X) doesn't depend on time. 



3 A possible example would be of the form X rn u^ n (x)uf n (y) , with appropriate decay of 

mSN 

\ m , regularity of u £ m of course. 
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Now recall the energy conservation for equation ([T]): 

(^Au^t) - V(x)u s (t),u s (t)) = (^Aug - V{x)u%,ul) => 



=> \\Vu^t)\\ L , < 2 V<-4A^ + v(,K. M g) + c.|| M g|| L ^ 

(75) 

Applying that to u|, it follows that there is a B > independent of A,e so 
that 

\\u\(t)\\H^B^±^ (76) 
for all (el. Now observe that (in the notations of lemma IB.1[) in general 

d Xi W e = Q[(d Xi + d Vi )p e ] 
(i.e. not only for pure states). Therefore (using equation (|166I0 

\\d Xi W E (t)\\ L 2 =e-%2n\\{d Xi +d yi )p e {t)\\ L 2 ^ 
^e-92nf\\(d Xi +d Vi )ul{x,t)lF x (y,t)\\ L 2d[x(\) ^ 
^ e-f-^Tr / (C\ + C,)dpL(X) s: B'e~ c 

for some B',c > 0. (We used the joint constraint on C\ : d^,(\) to interchange 
norms and derivatives with the integration). 
Now for the d ki derivatives, recall 

\\d kl W s (t)\\ L2 ^ \\d ki W%0)\\ L 2+27Tt sup \\d Xi W e (T)\\ L *; 

r£(0,t) 

see equation ([70)) for the derivation. The proof is complete. 

4 Proof of the main results 
4.1 Proof of Theorem QH 

The proof will be broken down to several steps. Denote be Wf the the SWT 
corresponding to the problem with the potential V\ , i.e. the solution to 



d t W( + (2vk ■ 3 X + e -^d x ■ d k ) W{- 

= °' (77) 



2 



Re 



i J e 2 ™ Sx -^ s *V x {S)W^{x + k - ?£)dS 



Wi(x,k,Q) = W$. 



23 



See equation ^ for the definition of V\. Naturally the corresponding Wigner 
function will be denoted by Wf (t) = <&~ 1 W( (t). Moreover denote by W e (t) 
the exact SWT satisfying equation ©, W e {t) = <&W e {t). Finally, recall the 
definition of p\, equation (1161) . and denote for further use 

n 

M,=V sup WltorKilWi&KrfWv*. (78) 
i= i *e[o,T] 

• Proof of equations dTTJ), (fl9|): 

There is a constant C, depending only on n,a x ,<7k, such that 

| \Wf - p\\\ L ? < CT{eHd 1 + y/iMi), (79) 
\\Wf - W e \\ L 2 CT{e^Mi + sM-l), (80) 
\\W £ (t) - W e (t)\\ L 2 <: Ce tD ^ n ^e\\W^\\ m , (81) 

and 

\\p\(t)-p{t)\\ L ^CT{e + ei)e TD ^\\W^\\ H ,. (82) 

(In the process of this last proof we will have to use that 9 > 1). Recall that 
D(l,n) was estimated in equation (|63[) . 

The conclusion then follows by estimating 

M x < ||Wf(t)|| ff x < e tD ^\\W^\\ m , (83) 

where theorem 12.21 was used for the last inequality - making use of Al(r=0). 

• Proof of equations (fl8|) . (f20]) : 

It suffices to combine equations (fTT)) . (flT?j) . theorem 12. 1[ and to use once 
again equation (j81) . to get 



\\W s (t)\\ L2 = \\W$\\ L i + y/iCe TD ^\\W§\\ H i (84) 

with |C| < C. 

Now we proceed to the proofs of the building blocks: 
Proof of equation ()79|) : 

One can Taylor expand W[ into a power series to recast equation ([77} as 

d t Wf + 2irk ■ d x Wf - ±d x Vi {x) ■ d k W{ = FWf (85) 

where 



F = -^fd x ■ d k - 



m=2 (m-i)mod2=l |A| = i 

|S| = m - i 



(86) 



4 because W e is an entire analytic function; see [3] for proof and more details. This is in 
contrast to the formal expansion often used for the Wigner transform. 
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and A,BG (NU {0})™ are multi-indices. See lemma POl for the full justification 
of this step (i.e. of interchanging the summation of the Taylor expansion and 
the dS integration). 

It is clear now that the SWT equation ([77)1 can be seen as a perturbation 
of the modified Liouville equation (|16p . We will use lemma B~2"l however before 
that some remarks are in order: 

Denote E\(t) the propagator of equation (|T6|) . i.e. p\{t) = £?i(t)Wo, and 
Ui(t) the propagator of equation ([77]). Wf(t) = Ui(t)W(j. The Liouville equa- 
tion (|16p can be solved with the method of characteristics for any e,r) > and 
its propagator has L 2 norm equal to 1, i.e. 

\\El(t)\\ L 2^ L 2 = 1. 

Moreover, it is clear that Ui(t) = $?7i(t)<i> -1 , where Ui(t) is te propagator 
for the corresponding Wigner equation, Ui(t)W^ = <I> _1 f/i(i)Wg . Following 
Theorem l2.11 U\(t) is an isometry in L 2 Mt E R. (Obviously a similar statement 
is true for U{t) : W$ i-> W £ (t)). 

Now by applying lemma R~2l it follows that, V0 ^ t ^ T, 



||/»f(t)-Wf(t)|| L2 <T sup HEi^-^FC/x^Woll^ < 

re[0,T] 

s;T sup ||F^f(r)|| i2 . 

r£[0,T] 



(87) 



Recalling equation we are called to estimate 

\\%f-d x -d k wf(T)+ 

+2 £ e + i «7?(-ir- 1 e ^£ (a) wwf(r)ii £3 < 

m=2 (m-l)mod2=l |A| = i 

|B| = m - I 

+2Eg E ^ E II^^^WfMlu, < 

m=1 (m-l)mod2=l \A\ = i 

|B| =to - i 

< W £ (t)||l2 + 

+2 g ^ E °* E lla "J^ llLX \\d^w[(r)\\ L2 . 

m=2 (m-l)mod2=l \A\ - I 

\B\ = to - / 

For a finite number of to's, namely to — 2 — 9 < —1, assumption Al(r=0) 
implies that 

\\d£ +B Vx(x)\\ L - <0(1). (89) 
As will be clear soon, these don't yield any intersting contribution. 



(88) 
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Now assume that m— 2 — 9 > — 1; it suffices to recall that by Al(r=0) 

||^ +B ^(z)lk~ < (27t)™m + D^^-^^r (ia^=i) , (90) 

where for brevity we denote 

/ n 2 

To estimate \\d£dj? W{ (r)|| L 2 we need to remember that \B\ ^ 1, therefore 
for each m, / appearing in the sum of equation (|88p there is a multi-index |6| = 1 
such that b ^ B. At this point one should recall the definition of Mi (equation 
(l78|) ). With that in mind, we have 

\\d£d%Wl{T)\\ L 2 =(2w) m \\X A K B e-^^ x \ 2 +^ K ^W{(X,K,T)\\ L 2 ^ 
^ (2Tr) m \\X A K B ~ b e-^^ x \ 2 +^ K \ 2 ^\ L ~\\K b W({X,K,T)\\ L 2 sC 

(91) 

< (2tt)(™" 1 ) ( ft ||^e-¥^|| ioo > ) ( ft ||^ fl - 6)- e-*^3|| ioo > ) Ml = 



d=l J \d=l 

A d (B-b) d 

(2^n(^) 2 (^) 2 



M 



(92) 



= Mi (|f)^^^ ( ) n AJ{B-b)—. 

d=l 

Finally, with the same arguments as above, observe that 

e -^\\d x ■ d k Wl(r)\\ L 2 < S|i27r £ \\X d e~^ x2 \\ L ~ M x = 

d=i 

Plugging equations (1T40|) . (I9T1) and $2$ in (JggJ we get 

\\FWt{r)\\ L 2 <^i^"M 1+ S 1+ 

m =[i+ei v 7 (93) 

n v (B-6) d 2 

2~i A\B\ 
(m—l)mod2=l \A\ = I 

\B\ =m-l 

Si stands for all the contributions from m — 2 — 9 < —1 (see equations (159")) , 
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(|90p ). Using our previous considerations it is easy to see that 

m* _. 

Ei = E &r E ^'ll^^WflUHI^Vx^IlL- < 

m=2 A,B 



(94) 



< CMi E & E ^' 'CU,b = O(VeMi). 

_ 2 

Since we already have a term of 0(y / eMi) - the one coming from ^-H&c • 
dkWi (t)||x,2 - we will not carry Ei explicitly in the sequel. 

Moreover, it is clear that (eventually) (27r) m M < ^^-(rl')- m= ^ =l T ( m ~^~ e ) 
so we don't treat explicitly the contributions from the first term of the rhs of 
l|90p for any m. 

Now, to proceed from equation (|93j) we have to understand better the inner 
sum. To that end, note that without loss of generality (specifically for m large 
enough) 

V V2ttA ' 

^ (95) 

r(^)<(2.)^2^y=^=:. 

The first statement is the Stirling approximation; if \A\ = it is clear that the 
contribution of the term A~? = 1 (this can be seen very clearly if one goes back 
to equation (|9"Tj) where it originates and check it. This point will always be 
understood as explained here, with no further explicit mention). For the proof 
of the second statement, see lemma l4~4l Now we have 



£g. {B ~ 2 b)d 

■ i ■ i 1 i d \^ ^0 d , - „ . 

E„l -r-m+l+l d=i p I m — 1 — 8 \ 

a x"k ATS! 1 V 2 I 

(m-l)rnod2=l \A\ = I 

\B\=m-l 



n f ; 2 - x /AMB-b) d \ 

_/ _-m + i+l d = l y/2*y/A d (B-b) d ^ ( m-X-j \ ^ 

2-^ a x u k A\B\ 1 V 2 7 ^ 

(m-X)mod2=\ \A\ = I 

\B\ = m — I 



\B\ =m-l 

(96) 



3+28 

, 24 

5> iin-l 



(27r) 4 ^/m(m-l) 



\ 



V ,_2i 2(/-m+l) ^ m! ^ 

"x "k A!B! 
1=0 \A\=l 

\B\ =m-l 



(2tt) 



3+2fl m— 1 

e 
2 
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where we used the observation that X Xd ^ \A^W E x d ano - the multinomial 

d=i y d=l 

expansion for the two last steps. 

Going back to equation (|93]) . we get 



(2tt)T— m =ri+e] V 



(97) 



, ,* M lD 2— ^ / e « 7r K + ^)\ 2 

— M £ h7 — a- 1 

(2tt)^^ m=2 \ 



.8 MiD2^^ A /n(l+o-|CT?) 
+7? 2 ; — - fcrr — ^— ' 



(27T) 



It is now clear that the constraints on r\ (recall that rf = r) -\ ^ , and r\ is 

specified in equation (|101| ) for this to work are 



Dl = -h^r 2-^<l 

»7H 



(98) 



^ = i^fe < 0(1) 

For example, by setting 

7T / 472 \ 

v > £ _ . max ^(4n - l)o* , _ _ CT 2 j (Q9) 

it follows that 

Dl < i, D 2 < -J-^. (100) 

In particular can be 

ry = e (f • max ((4n - l)o* , ^ - o*) + l) = 0(e). (101) 
Equation ({79")) follows by combining equations (1571) . (j!?T|) and (|101|) . 



Proof of equation (go]): It is obvious that \\W({t) -W E (t)\\ L 2 = \\$(W e (t)- 
Wf(t))\\L* < \\W%t)-W[(t)\\ L2 . 
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We will use again lemma l4~2l Recall the definitions of U(t), Ui(t). It is 
straightforward to compute that if 

LW{ = \Re\i J e 2mSx (l - e~ r >' s * )V (S)Wf (x , k - f)ds\ , (102) 

it follows that 

= f J Wf(X - S, K)(l - e- r '' s2 )V(S) [ e -^s-K _ e *ieS-K^ dS _ ( 1Q3 ) 

= 2nfWf(X - S,K) (l-e-v' s2 )V(S)S-K sm ^ s s K K) dS. 
Using lemma I4TT1 it follows that 

^2±n\\K l Wl{X,K)\\ LH ^ ) \\(\~-e^' s2 )V{S)S l \\ LHm \\ *&**p. \\ LaB . 
i=l 



(104) 

Obviously I "'"^s-'a^ I ^ ^' Moreover, recalling assumption (A3) we have 



" 1 

]T||i^f(X,K)|| L2(R2n) < Mi — . (105) 

Finally, recalling equation (fT2|) of Al(r=0) 

||(1 _ e-v' s2 )S t V{S)\\ LHm < 0(r/§) (106) 

without loss of generality (i.e. 6* < 2 without loss of generality). 

Therefore summarizing equations (|104[) . (|105[) . (|106[) we finally get 

\\LWl\\ L , { ^ n) ^M.Crfi. (107) 

Now by applying lemma l4~2l we get 



t 

W({t) - W £ {t) = - J U(t- t)LUi(t)W§(1t =>■ 

T = 

=► ||Wf(t)- W*(t)|| L a < 

sup \\U(t)\\ L 2^ L 2 sup \\LWf(t)\\ L2 < 

re[0,T] re[0,T] 

s: t (y ^ + ??') C" Mi 
for some 0(1) constant C . 



(108) 
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Proof of equation (|8ip: Assumption Al(r=0) allows us to use theorem 12.21 
for m = 1: there is an O(l) constant D(l,n) so that 

^e tD ^\\W£\\ Hl . (109) 

Now the conclusion follows by application of lemma HOI - and more precisely 
its corollary described in equation (I16ip . 

Proof of equation (|82[): For this we will introduce an auxiliary object. Denote 
P2{x,k,t) as the solution of 

d t p 2 + 27rfc • d k p2 - ^d x V ■ d k p 2 = 0, 

(110) 

P2 (t = 0) = W§, 

First of all observe that, using again equation (I161|) . there is a constant C, 
depending only on <r x , a^, n, M , such that 

\\ P2 (t) ~ P (t)\\ L 2 = \\wz-ws\\ L i ^c^\m\\ H i. (in) 

We will proceed in a way analogous to what we did earlier to show that p 
and pi are close. 

By substracting equation (|110j) from (fTBj) it follows that h = pi — p 2 satisfies 
d t h + 2irk ■ d k h - ^d x Vi ■ d k h = 

= ^x.k^A^ I (1 - e-v's 2 )V(S)p 2 {X - S, K)S ■ KdS], (112) 
h(t = 0) = 0. 

By using again lemma [4T2l - it is clear that Liouville equations with W 2 '°° 
potentials are well posed in L 2 - it follows that, for t € [0,T], 

\\h(t)\\v < T||(l - e-Vs'jIfl'lVCfiOIUiMi < TMiCCe + el). (113) 
Recalling equations (|106l) and (1531 . equation ([52]) follows. 

The proof of theorem ll.2l is complete. 

4.2 Proof of Theorem Q 

The proof is to a large extent analogous to that of theorem ll.2l fusing the H~ m 
continuity of U (t) , theorem 13.31 instead of the L 2 continuity - theorem 12.11 - 
which we used before). Therefore, we will not discuss in detail points that were 
treated previously. 

The main idea is that 

\\W £ (t) - p(t)\\ H -r < \\W £ (t) - W S {t)\\ H -r + 
+ || Wf - WZ\\ H - r + \\Wf - pl\\ H -r + \\ P \(t) - P (t)\\ H -r. 
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It will be shown that there is a constant C such that 

Wt - Pl\\u-r < CT\\W^\\ H -,-re TD ^ n \ei + y/e), (114) 

\\Wl-W e \\ H - r <CT(e* +e)\\W^\\ H -i-re TD( - r+1 - n \ (115) 
\\W e (t) ~ W e (t)\\ H -r < C\\W s \\^-re TD ^- l ^yfe, (116) 

and 

M{t) - P{t)\\n-r < 
< {e + ei)CTe T W r+1 ^+ D ^\\W^\\ H -i-r+ (117) 

+c^ TD{r > n) \m\\ m ~ r . 

Now we proceed to the proofs of the building blocks: 
Proof of equation (|1 14|) : 

Recall that E\(t) is the propagator of equation (TloT) . i.e. p\{t) — E\{t)WQ, 
and Ui{t) the propagator of equation ([77)1. (t) = JJi^WJ. 

It follows from theorem 12. 31 (applying the second sufficient condition), using 
the assumption Al(r), that 

\\Mt)\\H^Hr ^e w ^ n \ 
With respect to the sufficient condition, observe that 

sup H^lli- < ||Fi(5)|^r +2 || il < oo. 

|a|^r+2 

Therefore, it follows by duality (just as in lemma I3~3l that 

\\Ex{t)\\H-r^H-r^e tD ^\ 

(Recall that D(r,n) was defined in equation (|62l) . and can be estimated using 
M(r)). 

Moreover, it follows from lemma |3~31 that Ui(t), U(t) are bounded in H r -> 
H- r . 

Now by applying lemma l4~2l it follows that, VO ^ t ^ T, 



M(t)-Wl(t)\\ H -r S U P \\E X {t-T)FU X {T)W§\\ H -r < 

re[o,T] 

_ (118) 

^ Te TD(r,n) gup 1 1 .F Wf (r) 1 1 ff -r . 
re[0,T] 



31 



Recalling equation (|5o]l . we are called to estimate 
\\^d x -d k W({r)+ 

E i<-™+V«(-ir-< e ^^W^fMlli 



(IFF 

m=2 (m-l)mod2=l 



A = i 
|B| = m - i 



2 

+2E^ ^ ^ E ll^§^WW?(r)|U-, < 

m=2 (m-l)mod2=l \A\ = I 

\B\ =171 — / 

^ s -£\\d x -d k Wf(r)\\ H -r + 

+2 E E E ii^ + ^ r )iu- pA^ iS(r) || H _ r , 

m=2 (m-!)rood2=l |A| = I 

\B\ = m- I 

_ (119) 
Now we proceed to estimate \\d x d B Wl (t)\\h-t-: 

\\d£d?w;(r)\\ H -r = 

= (2Try n \\X A K B e-^^ x \ 2 +^ K ^Wf(X,K,T)\\^ H - r < 
< (2n) m \\X A K B - b e~^^ x ^ + ^ K ^^\ L ^\\K b Wf(X,K,T)\\^ H - r < 

<(27r)(m -i) (njix^e-^^iuoc) 

d=l / 

tD(r+l, 



= HWJHff-!-^" 3 ^) (%) m * 1 v?vJ m - l - 1) n V (5-6)7 

d=l 

(120) 

Finally, with the same arguments as above, observe that 

E -£\\d x - d k W?(T)\\ H -r < 

< ^2tt £ ll^e-^^^H^ ||WJ|| H -i-re tD (H-i.») = ( i2l) 
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Plugging equations ([IT]). flTJO} and (fT2H in (fTT^|) we get 

\\FWi(r)\\ H -r < ^ ^||W^ ( f|U- 1 -. e T ^+ 1 ^ + S 1+ 

m =ri+e] v ' (122) 

n AZ(B-b)^- 

{rn-l)mod2=\ \A\ = I 

\B\ =m - I 

Si stands for all the contributions from m < + 1, and behaves simlarly to 
v ^||W £ || // -i-.e T£ '( r+1 '™) (see e.g. equation dH ). 

By working in complete analogy to the proof of equation (f79|) . (|114j) follows. 

Proof of equation (|115|) : It is obvious that 

\\Wm ~ W B {t)\\ H -r < \\Wm - W S {t)\\ H -r. 

We will use again lemma [4T2l to estimate 1 1 WJf (t) — I 7 ^ (£)!!#->-. It is straight- 
forward to compute that if 

LW{ = lRe\i f e 2mSx (l - e- r >' s2 )V(S)Wf(x,k - ^)ds] , (123) 

it follows that 

\\LWf\\ H -r = \\?[LWfi\\j: H -r < 



<2J2Tr\\K i Wl(X,K,t)\\r H -r \\(l-e-"' s2 )V(S)S t \\ Ll \\ ain E % s K K) 



2n7r||W £ ||H-i-'-e T - D(r+1 ' n) M 3 . 

(124) 

Recalling equation ^ of Al(r=0), ||(1 - e~ r '' s2 )V(S)S i \\ L i = 0(e + el); 
equation ()115|) now follows. 

Proof of equation (jll6[> : It follows readily by application of lemma IA. II for 
p = 2, m = — r, s = 2. 

Proof of equation (|11T[) : Recall the definition of p2(x,k,i), equation (jl 10|) . 

First of all observe that, using again lemma IA.11 there is a constant C, 
depending only on <j x , <Jk,n, M 0} such that 

\\p 2 (t)-p(t)\\ H -r < e TD ^\\W$ -WZ\\ H - r < CV~ee TD ^\m\\ m - r . (125) 
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By substracting equation (|110j) from (fT()|) it follows that h = p\ — p2 satisfies 
d t h + 2nk ■ d k h - j^d x Vi ■ d k h = 

= J^K-t^JO- ~ ^ V ' S2 )V(S)MX - S, K)S ■ KdS], (126) 
h(t = 0) = 0. 

By using again lemmata 14.21 14. fl it follows that, for t £ [0, T], 

\\h(t)\\ H - r <CTe™('-.»)||(l-e-'''^)|5|y(5)|U 1 ||p 2 (t)|| H - 1 -. < 

(127) 

«S (e + e§ )c'Te T ^ r+1 ^+ D ^ \ \ W§ \ \ H -x-r . 

for some constant C depending only on n, a x , a kl Mo. (We also used theorem 
231 to estimate the rhs of \\h(t)\\ H -r < \\pi(t)\\ H - r + \\p 2 (t)\\ H - r ). 
Using the triangle inequality, equation (|117|) follows. 



The proof of theorem 11.31 is complete. 
4.3 Proof of Theorem 11.41 

The proof is very similar to that of theorem ll.3l We will only discuss explicitly 
the points where differences arise. 
Of course the main idea is that 

\\W e (t) - P (t)\\ Hr ^ \\W"(t) - W S (t)\\ H r + 
+ - W^\\ H r + \\Wf ~ pi\\ H r + M(t) - P (t)\\ H r. 

It will be shown that there is a constant C such that 

\\W[ - pl\\ H r ^ CT\\W^ + re TD ^ +1 ^[ei + Vi), (128) 

\\W[ - W e \\ Hr ^ CT(e + e S -^)\\W^\\ H i + ,e T[D(r+1 ' n)+D ^ n) \ (129) 
\\W*(t) - W%t)\\ H r < C\m\\ Hl+r e TD ^ +1 ^V^, (130) 

and 

M(t) ~ P<t)\\HT < C^e TD ^)\\Wt\\ m+r + 

+CTe T[U(r,») + D(r+l,n)] 1 1 W * 1 1 Jf „ +1 ( e + e V ) . 



(131) 



Now we proceed to the proofs of the building blocks: 
Proof of equation (|128|) : It is entirely analogous with the proof of equation 
(|114p . in the proof of theorem ll.3l (section l4.2l) - it is not necessary to replicate 
it here. 
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Proof of equation (|129|) : It is obvious that 

\\Wf(t) - W S (t)\\ H r < \\Wffl - W S {t)\\ H r. 

We will use again lemma 03] to estimate \\Wf (f) — W £ (t)\\H^ ■ Using lemma B~T1 
it follows that, if L is defined as in equation (|123[) . then 

1 1 LWf | \ H r = \\T[LW{]\\r H r <i 
<C|| \K\Wf(X,K,t)\\„r 11(1 -e-^ 2 )(\S\^ + \S\)V(S)\\ L , \\ sin ££ K K) \\ L « 
< C\\W§\\ H . +r e TD ^n) _ e -^)(| 5 |H-i + \ S \)V (S)\\ L ,. 

(132) 

Now recall assumption Al(r) to see that 

||(l-e-^ 2 )| S |r+iy (5) || il(Kn3<C - (£ + e ^) (133) 

Therefore 

\\LW[\\ H r < C(e + e i i r )||Wo e || ff i + re TI3 ( r + 1 '") (134) 
Equation (II 20[) follows by applying lemma [ 



Proof of equation (|130p : It follows readily by application of lemma IA.1I for 
p = 2, m = — r, s = 2. 

Proof of equation (|131[) : Recall the definition of p 2 (x, k,t), equation (jl 10|) . 

First of all observe that, using again lemma IA.11 there is a constant C, 
depending only on a x , Ukn, M , such that 



\P2 



(t) - p(t)\\ H r < e TD ™\\W§ ~ WZ\\ H r < CV~ee TD ^\\W^\\ m+r . (135) 



Applying lemma |4~T1 to equation (|126l) once more (and using again Al(r) ), 
it follows that for te [0,T] 

||M*)lk' <C\\(l-e^' s \\S\ + \S\^)V(S)\\ L1 \\p 2 (t)\\H^ < 

(136) 

< Ce TD ^ +1 ^\\W^\\ Hr+ i{E + e a ^r). 
Using the triangle inequality and lemma [4721 equation (1131j) follows. 

The proof of theorem ll.4l is complete. 

4.4 Auxiliary lemmata 
Proof of lemma 11.11 
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First of all 

J\V(S)\\S\^dS^ J \V(S)\dS+ J \V(S)\\S\^dS ^ 

\S\<1 \S\>1 

oo 

< Mo + D J s r+2+n-l-n-l-6 > ds < ^ 
s=l 

because r — < — 1. This, together with equation (fT3"|). implies equation ((HJ. 
Now denote for convenience 

v' = V+^- (137) 

One observes that 



m+ B Vi(x)\\ L oo < (27r) m / ins) (n ^ A+fl)i ) e""'l s l 2 |^ < 

<(2^r/|T/(5)||5r" e -"'l s l 2 ^< (i3 8 ) 
(2tt)" 1 ( / \V(S)\\S\ m dS + D J r (™-2)-e e ->n'r 2 dr j _ 

V|S|<1 r=l / 

For a finite number of m's, m — 2 — 6 < — 1, in which case, using equation (|1 38[) 
and equation (|8|), for |A| + \B\ ^ m we have 

ll# +B Vi(*)IU~ < (27T)™ M + D / r( m - 2 '- e dr ] = 0(1) (139) 



r=l 

Now assume that m — 2 — 9 > — lp then equation (|138[) leads to 



||^+ B Fi(a;)|| L oc < (2^) m M + Z?(27r) m / r (m-2)-f e - n 'S ' d5 = 

= (27r) m M + Dteft-tfy^^r {^^) . 

This concludes the proof of equations (|10l) and (fTTI) . 



(140) 



5 Wc don't have to worry about m — 2 — 9 = —1 because we have restricted 9 £ N. 
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Now observe that 

||(l- e -Vs 2 )| 5r+ iy ( 5)|| il(Rn)< 



< / {l-e-v' s2 )\S\ r+1 \V{S)\dS + D J \S\ r + 



|S|>1 



\S\r- 



8=1 



< r/M + L> 



fs r - 1 - e rfa 2 d8+ J a 1 - 1 - 9 da 

1 a=vi 



=: M a . 



By selecting A — rj' 2,we get 



U,(r l = ,/UM-/)(r ? ^^^; e - 1 + ^ 



By recalling our assumption that 6* > r + 1 it is easy to see that 
M 3 (r) s? C(£^ +e) = o(l). 



(141) 



(142) 



(143) 



This concludes the proof of lemma 11.11 

Lemma 4.1 (A version of the Young inequality) Let f G L 1 (R"). Then, 
ifgeL 2 (R 2n ), 



II / f(s)g(x- SjfcJdsHiJpRJn) < ||/||il(Kn)||fl , |U = (ltt=™); (I 44 ) 

«er m (i 2 "), 

II j f{s)g(x - S,k)ds\\jr H - m{S 2 n) < ||/||Ll(R»)||S'||^H-m( M 2n); (145) 

and i/je if m (R 2n ), /(s)(l + |s| m ) G L 1 (M n ), 

/m 
/(s)fif(a; - s, fc)^!!^™^) < CllffH^Hm^n) ||/(s)|s| i |ki( R "). (146) 



2=0 
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Proof of equation (|144|) : 

II / f( s )g( x - s, k)ds\\ 2 L2{R2n) = J f(s)g(x - s, k)f(y)g(x ~ y, k)dxdk 
I j g(x - s, k)g(x - y, k)dxdkf(s)f{y)dsdy 1 1/| || 1(Kn) \\g\ |£ 2(K 2„). 

Proof of equation ({145)1 : Denote B = {(f>\ \\(j)\\^H^ < 1}. Then 

II / f(s)g(x - s,h)ds\\ H - m( yp n ) = sup J f(s)g(x - s,k)ds(j>(x, k)dxdk < 



>(R") SU P J g(x - s,k)<p(x,k)dxdk < 1 1/| jr,i( R ") I Iffl I j?- 



Proof of equation (|146|) : Recall that 

II I f( s )g(x - s,k)ds\\jr Hm{R 2 n) = 
= E \\(2irx) A (2nk) B J f(s)g(x-s,k)ds\\ L 2 (M 2 n) , 

\A\ + \B\^m 

therefore it suffices to work for one choice of A, B with \A\ + \B\ ^ m: 

11/ f(s)g(x - s,k)dsx A k B \\ L 2 (R 2r l) = 
= \\I f (s)g(x - s, k)ds(s + (x - s)) A k B \\ L 2 {W 2 n) < 

< E (^H/^ls^-'lUH^IIfl^^^^lU^M-) < 

(=0 

^\\g\\H^)£ti)\\f(s)\s\\ A -n LHRn y 

1=0 

Lemma 4.2 (An observation on perturbations) Consider two equations, 

vir ' ( i47 ) 

v(t = 0) = w , 

and 

u t +Tu = Lu, 
u(t = 0) = wo- 

Denote E and U their propagators respectively, i.e. v(t) — E(t)u>o, u(t) = 
U(t)wQ. Then the difference of the solutions can be cast as 

t t 
u{t) - v{t) = - J E(t- T)LU{T)w dT = J U(t - t )LE(t )u>o dr (149) 

T = T=0 

Moreover, if 

v T ir m ( i5 °) 

r(t = 0)=w , 



then 

t 



r(t) - v(t) = - / E(t- T)f(r)dr. (151) 



Remark: The considerations here are formal; questions of well-posedness of 
the equations etc should be considered each time the lemma is applied. 
Proof of equation (|149|) : Set 

ip(t) = U(-t)(u(t)-v(t)). (152) 

Then 

U{i)ip(t) = u(t) - v(t) U(t)ip(t) + U{t)ip{t) = u{t) - v(t) =>> 

=> -(T + L)U{t)ip(t) + U{t)ip{t) =-Tu[t)- Lu(t)+Tv(t) 

— (T + L)(u(t) — v(tj) + U (t)ip(t) = -Tu(t) — Lu(t) + Tv(t) => 

t 

=> tjj(t) = -U(-t)Lv(t) U{t)i)(t) = -U(t) f U(-r)LE(T)w dr 

t 

=> u(t) - v(t) =- J U(t- t)LE{t)w <1t. 

T = Q 

(153) 

The other equality of equation (I149P follows essentially in the same way. The 
proof is complete. 

t 

Proof of equation (fT5T|) : Set r'(t) = v(t) - f E(t - r)/(r)dr; by straight- 

T = 

forward computation it follows that dtr 1 + Tr' = f(t). 
The proof is complete. 

Lemma 4.3 (Interchanging summation and integration) We will show that 
| e 2 II ^f^Y 1 (S)^( I + ^,^f) ( |S = 

9* 9 ~ +°° m n ■ 

= /e Wl - ( ' + ^" )s V(S) £ ^[Y,Sd(i<T 2 x d Xd -d kd )rWl(x > k)dS = 

m=0 d=l 

+ 00 9 9. ™ ■ 

= E fe^ s *-^+^*)s v^)-^^ S d {i°ld Xd - d kd )TWt{x,k)dS, 

m=0 d=l 

i.e. that interchanging summation and integration is justified. 
Proof: Denote 

g(S) = e 2°* s Wt{x+^-,k-—), 
9m (S) = e-^°l)S 2 g -^[J2s d (i*ld Xd - d kd )] m W[(x,k). 

m=0 ' d=l 

We know that Egm('S') = d(^) pointwise (even uniformly on compacts), 

m 

since it is the Taylor expansion of an entire analytic function. 
Observe that since the potential always satisfies Al(r=0), 

V(S) £ L 1 ; (154) 
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therefore, it suffices to show that 

£>m(S)| <<7, 
m 

the result then follows by dominated convergence. Indeed, 

\9m(S)\ < 

+00 n a 

m=0 |A|=m d=l 

„ +00 n Ad 

<C e -('7+¥^)'5 2 ^ e m V- Cm\ oA Tt (A 

ra—\ 

Here observe that 



=0 \A\=m d=ll=0 



(155) 



(ia* x d Xd Y(-d kd ) A «- l Wt(x,k) 
< m aS {f,CT2i}(2 7 r)™/X^ d -'e-¥(^l^l 2 +^l^l 2 )|iy 1 (X,X)|dX ( iX 

Moreover, 

\\X l d e-^\ x \ 2 \\ 2 L2 = (Je-^^dx)^ 1 Jx 2l e- S ™> 2 dx = 



HE; 

(7T£0-J) 2 



= v " — ■ 



and similarly 



and therefore 



\\K^- l e-^l\ K f\\l 2 = {ealf 



^ 2(A d -Q-l ^ 



L 2 ~ \ CU kJ ' „ 2(A d -Q-l ! 



1-21 l-2(A d -Q 



< max{l, a 2l }(2Tr) m (ea x cr k ) 2 ( e7r ) 2 cr^. 2 a 

^r(M)r(2^ 11^11,2 < 

1-21 l-2(A d -Q 

^ Cmax{l , a x }max{l , <7 X 2 }max{f,<r fe 2 } 

(2.)- («r) ^ /r(2i_i)r(^^) < 



< CB m e—JT 
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where C, B < oo are independent of m. Plugging this back into equation (|155[) , 
there are (possibly new) constants C, B < oo such that 



| 5m (5)KCe-(^T^)s 2 Y. B 

n A d 



m e 

■m ! 



E G)^nE(1^Jrt^)rp^Hk 



1-4.1 = 



2=1 /=0 



-f oo 



n A d ./r 



m=0 |A|=m d=li d =0 



Zd!(A d -Z d )! 



< 



< C'e-^+^-Ds 2 + f B m ef E ^ ft E ^SSF^ ' 
At this point, using the Cauchy-Schwarz inequality one observes that 



(156) 



fr ^ VCi-l)'(Ad-id-l)! _ fy / ^ A,! 

11 ^ id!(A d -Z d )! - lM/ 2^ 



Ad-1 



d=i V id=o 



idl(A d -J d )! V^TA/ E ; d (A d -z d ) 



< n 

<Z=1 



A d 



E 



A rf ! 



A d !(A d -l) \ l d \{A d -l d )\ ~ V Ad! ' 

V td— 



Finally, observe that 
so altogether (possibly for new constants C, B) 



f7T 



|ffm(5)| < C E B m £T(ri + fal)-^ E T^y < 

m=0 |A|=m ' 



e ^(i + f^)-^ E 

m=0 |A|=m V 



+ ■00 



B™eV(r,+ fal)-^l E 



|A|=m 



mW |A| = n 



(157) 



+00 



< C E 5 m e~(»? + faD-^Vrf' 

m=0 



+ °° (H'\» 

m— 



(For the last step it is important to recall that finally rj is fixed r\ = Bq£ for 
an appropriate constant Bq, see equation (jlOip 1 ). 
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This justifies the interchanging of the integral and series. To proceed to 
equation (|86|) then we carry out the (inverse) Fourier transforms in 



^ )s ~v(s)^[Y. s d (i<r 2 x d Xd - d kd )] m w[(x,k)ds = 



= T. 



v(s)^[E s d (i* 2 x d Xd - d kd )rwax,k) 



The proof is complete. 
Lemma 4.4 (An observation on the Gamma function) Form large enough 

' m — 1 — i 



,„ . 1 _ 3 + 28 

< {2tt)±2 — 



2 m m(m- 1) 



,we re 



Proof: The Gamma function satisfies the duplication formula, 

T{z)T(z +\) = 2^- 2z V2^T(2z). 

Keeping in mind that T(m - 1 - 9) < T(r7i - 1) = (m - 2)! = m( ^!_ 1) 
working for m large enough), and setting z — T "~ 9 1 ~ 6> it follows that 

r(s=£=£)r(2i=£) = 2i- m + 1+e V2^T(m -1-6) 

(r(Hi^«)) 2 < T( r -^^)T{^) = 2i- m+e V2^r(m - 1 - 9) < 

< 2i+ e 72T 9m "' ! 

v 2 m m(m — 1) 

The proof is complete. 

Lemma 4.5 (Sobolev norms of delta functions) Let z$ G R 2 ™. Then 

2m >n => 5{z - z ) e # _m (R 2 ") 
Proof: Take / e H m (M 2n ). Then 

\(6(z-z )J(z))\^J\f(x,k)\dxdk = 

= J \f(x,k)\(l + \x\) m {l + \k\) m (l + \x\)~ m (l + \k\)- m dk ^ 
y/j |/(fc)| 2 (l + |a;|) 2m (l + \k\) 2m dkyjj (1 + \k\)- 2m dkf (1 + \x\)- 2m dk < 

<C||/|| H ™ (l + + f p^-^dp J 



The proof is complete. 
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A Two useful properties of the SWT 

Lemma A.l Recall the definition of the spaces X m ' p , equation |ff7| j. Then, for 
any m G K, p G [1, oo], s > 1, the following estimate holds 



\W - W\\ xm „ < (1 + -max{^, a^}) e * ||W|| .^a,, 



(158) 



Remark: The statement and the proof are optimized for a x ,ak = 0(1). When 
that is not the case, fundamentally the same approach can still be made to work. 
We don't treat that here explicitly. 

Proof: For p G [1, oo), we have 



W-$W\\% m , r = \\W(X,K)(l-e-^^ x >)\(X,K)\ 



m\\P _ 
I LP — 



J (1 - e~^< x +< K ")\(X, K)\ m W{X,K) 



J (1-e-Tl^ 2 +< K ^)\(X,K)\ m W(X,K) 



I (1 



V £ 



e 2 



i K i)\(X,K)\ m W(X, K) 



dXdK = 

dXdK+ ( 159 ) 
dXdK < 



^(imax{alal}) p s^ J \(X, K)\ m+ ~W(X, K) 

X 2 +K 2 <\ 

+ / \(X,K)\ m W(X 1 K) P 



dXdK- 
dXdK < 



^^max{alal}) p ei\\W\f xm+ ^ + 



\(X,K)\ m W(X,K) 



^X 2 +K 2 >-^ 



dXdK. 



for any s > 1. 

For ^2 it suffices to observe that 

W\\ p m+2 , p = / \{X,K)\^+i)p\W(X,K)\PdXdK > 



J \(X,K)\ m P\W(X,K)\PdXdK = e-f/ 2 . 



(160) 



For p = oo it follows along similar lines (the only difference is technical, we 
can't raise to the p — oo). 
The proof is complete. 
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It might be instructive to write down explicitly some corollaries of lemma 

ED 

\\W-W\\ L * ^(l + ^max{<T 2 x ,*t})ei \\W\\ H i, (161) 

\\W- W\\v < (1 + ^max{a 2 x ,a 2 }) e \\W\\ H 2, (162) 

\\W-W\\ H -^(l + ^max{a 2 x ,*l})ei \\W\\ L ,, (163) 

\\W- WWh-! < (l + ^maxK,^}) e \\W\\ H i. (164) 

Lemma A. 2 (Critical smoothing) If a x a k = 1, then 
W e (x,k) = (Op Weyl (W s )(t>l tk ,^, k ), 

where 




If a x <Jk = 1 we will say that we have critical smoothing, and the SWT 
coincides with the Husimi transform. 

Moreover, if Opweyi (W E ) is a nonnegative definite operator, and cr x o~k ^ 1, 
then W £ ^ 0. 

Proof: Recall that 

(L(x, k), W s [f, g])=f L(x, k)e 2 ™ k yj(x + f)g(x - f)dydxdk = 

= J L(^,k)e 2 " k ^ a dkg(X)dXj(Y)dY = 

= / L(^,ek)e 2 ^ x - Y Ukg(X)dXj(Y)dY = {Op Weyl (L)gJ). 
Moreover, it is straightforward to see that, as long as a x ak = 1, 



F XoM (x,k) = 




= W s [<jf XoM ]{x,k). 



Therefore 

W £ (x ,k ) = (W £ ,F XoM ) = (Op Weyl {W s )4>1 m^1 m)- 

For the second part of the theorem, it is obvious that if a x a k > 1, W £ can be 
cast as a Gaussian molification of a Husimi transform. But the nonnegativity of 
the operator passes on to the Husimi transform, and any Gaussian smoothing 
preserves it. 
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B On the initial data of the Wigner equation 

In the case of a pure state W £ [u e ](x, k) = J e~ 2 ^ lky u £ {x + ^f)tiF(x - ^f)dy the 
following identity holds: 



J W £ {x,k)dxdk = ||w e ||| 2 . 



More generally, in the case of a mixed state W e {x, k) = j e~ 27Tlky p(x + x — ^f)dy 
we can recover the trace of the operator with Weyl symbol W - equivalently, 
with kernel p: 

W £ (x,k)dxdk = J p{x,x)dx = tr (OpwW(x,k)) . 

This means that the conservation of J W £ (x,k)dxdk in time evolution re- 
flects a significant physical fact; it emerges therefore as a reasonable normaliza- 
tion to set 



W £ (x,k)dxdk = 1 (165) 

for all e > 0. This is the normalization consistent with the standard scaling of 
the Schrodinger equation, and the concept of a Wigner measure |12) , 



W £ 



W°, J W £ (x,k)dxdk = J W°(x,k)dxdk = l. 



At the same time, the L 2 norm in also preserved (see theorem 12. and this 
is also of physical significance, 

\\W £ [u £ ](x,k)\\ L2 = e -f||u e ||| 2 . 

In the case of mixed states it corresponds to the Hilber-Schmidt norm, 

\\W £ (x,k)\\ L 2 =e-*\\p(x,y)\\ L3 =e-S\\OpwW(x,k)\\ H s. 

Since we deal with linear homogeneous equations, switching between the 

two is completely painless. The results here are formulated without any one 

normalization being fixed, that is why the assumptions involve quantities of the 

form H ty °H"'" 
tOTm \\W§\\ L 2 ■ 



It is instructive to look at pure states and how they scale: 

Lemma B.l For the Wigner transform of a pure state, W £ = W £ [u £ ], the 
following holds: 

\\d Xz W £ \\ L 2 <2 £ -f \\d Xz u £ \\ L 2 \\u £ \\ L 2, 

\\d kt W £ \\ L 2 ^47T£-f \\fu £ \\ L 2 \\U £ \\ L 2. 
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Proof: Denote by Q the mapping 

Q : f(x, y)^ J e^fix + f , x - f)dy. 

For example, the Wigner transform of a pure state can be recast as 

W £ [u £ ] = Q[u £ {x)u £ (y)\. 

It is easy to check that 

IK2/II^=^II/IU- (166) 

Moreover, it is easy to check that 

d Xt W'[u e ] = Q [(d Xi + d Vi )u £ (x)u £ (y)} , 

d kt W £ [u £ ] = Q [-2m^±u £ {x)u e 
The proof is complete. 



The importance of the following lemma is that it assures us that for coherent 
states the convergence to the Wigner measure is in fact in H~ r norm, not just in 
weak-* sense. This is in particular useful when working with strong semiclassical 
limits for pure states. 

Lemma B.2 (Wigner measure for coherent states) Assume Uq(x) is a co- 
herent state, 

u s (x) = e-faf^Je 2 "^ 1 , (167) 
v e 

with ae5. 

Then, for $ £ S, and for all r G N, 2r > n, there is a constant C depending 
on a(x), r such that 

\(W £ [u £ ] - 5{x ,ko),4>)\ < C\\<p\\ H r+2 +e ^) . 

Remark: Similar results can be shown for WKB initial data, i.e. u e {x) — 
A{x)e 2 i ±s( - X ) . 

Proof: It suffices to work for xq = ko = 0. Indeed, 

(W £ [u £ },<f>) = e-3 / e - 2 ^a(^3-)a(^ 7 ^)(/)(x, k)dydxdk = 

= J e- 2mk y a (x + ^-)a{x - ^-)<p(^x, k)dydxdk = 
= f e- 2niky a(x + ^)a(x - ^)<f>(0, k)dydxdk + 

J 2 Z 



+ f e- 2mk y a {x + ^-)a{x-^)R 1 dydxdk, 
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where of course R\ is the remainder of the Taylor expansion. We will collect 
and estimate all the errors in the end; in the meantime we have 



/ e- 27Tlk y a (x + 4 E )a(x - ^-)(j>(0, k)dydxdk = 
= J e~ 2vikv 4>{Q, k)dka(x + &*-)a(x - ^-)dydx = 
= Jfa(Q,y)a(x + &)a(x-&L)dydx = J fa(0,y)[a(x) + E{\[a{x) + E[]dydx = 
= 0(0, 0) + / 02(0, y)[a(x)E 1 + a{x)E' 1 + E^dydx, 



where of course E\, E[ are the Taylor remainders from Taylor expanding a(x± 



^) around x, and <p 2 (x,y) = J e 27viky <p(x, k)dk. 

Now, about the errors: \R\\ < |a;|\/fF||V</>||z,°o, an d therefore 

|/iK V^\\V<f>\\L~ \ J e- 2mk y a (x + ^t)a(x - ^)\x\dydxdk\ < 

< V^HV^Iioo \jS(y- 0)a(x + ^±)a(x - y^l)\x\dydx\ < v^I|V0||l~ / a 2 (x)\x\dx. 
Moreover, < ^\y\y/s\\Va\\ Lao . Therefore, 

\JM0,y)a(x)E 1 dydx\ ^ v f^ \jfo{0,y)a(x)\y\d y dx\ = 

= ^N V 2 ° I I*" / 102(0, y)\ \y\dyf\a(x)\dx. 



Finally, for the EiE[ term, one observes that 



/ M0,y)E 1 E' 1 dydx = 
= J 02(0, y)[a(x + ^)a(x - 4 E ) - a(x)E 1 - a(x)E[ - a 2 (x)]dydx 



and since we can select a M > 1 such that 




02(0, y)[a{x)E 1 + a{x)E[ + a 2 (x)]dydx < \fs / $2(0, y)dy, 




it follows that 



\jM0,y)E 1 E' 1 dydx\^\ J M0,y)E 1 E' 1 dydx\ + 



/ \M0,y)\dy + \ J M0, y )a(x+^)a(x-^)dydx 




Now we have 



I / M0,y)E 1 E' 1 dydx\ ^ ^\\\7a\\ 2 Loa J dx J 2 (O , y)\y\ 2 dy , 
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h~\ J MO,y)a(x + ^)a(x-^)dydx\^ 

\x\>M, j/EM" 

<| J MO,y)a(x + ^)a(x-^)dydx\ + 

|x|>M,|j/|>AT 

+ | / MO,y)a(x+^)a(x- ^-)dydx\ < 

|x|>M, |y|<M 

<| / MO,y)(fa(x+^)a(x-^)dx)dy\+ 

\y\>M 

+ | / MO,y)a(x+^)a{x-^)dydx\. 

\x\>M, \y\<M 

At this point observe thal[^| 

|z| > M, \y\<M =► < |ar ± ^| < 2|z|; 

now making use of the fact that a e 5 it follows that, as long as |x| > M, \ y\ < 
M, there are C 2 > 0, si > n such that 

\a(x + — - )a(x — )| < C 2: 



2 ' y 2 71 " (l + M) Sl 
so that finally 

h<\Hx)\\\, J \M0,y)\dy + C 2 J \fa(0,y)\dy J 

\y\>M \y\<M \x\>M 

To summarize (and keeping all the M dependencies explicit - that will be 
important later), 

\(W s [u s o ]-S(x-0,k-0),<t>)\ < 
^\\a(x)\\l 2 J \M0,y)\dy + C 2 J |&(0, y)\dy M n ~^ + 

\y\>M \y\<M 

+Ce\\Va\\l ! „M"JMO,y)\y\ 2 dy+ Vg||V 2 a|k °° / 1^(0^)1 \y\dy f \a(x)\dx. 

(168) 

Obviously 

/ 102(0, y)\dy < / \<j>(x,y)\dxdy, 

\y\<M 

JM0,y)\y\ 2 dy^f\j(x,y)\ \y\ 2 dxdy, ^ 
J\M0,y)\ \y\dy^f\$(x,y)\ \y\dxdy. 



by restricting, without loss of generality, e € (0, 1) 
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Now it suffices to see that 

J\$(x,y)\(\l+y\) 2 dxdy = 

= / \$(x, y)\(\l + l/lPp+^^ydl + x\Y(\l + y\Ydxdy < (170) 

< C\\(t)\\H^A\ (\l + X \Y(\ 1 + y\)A \L^ = C'\\(t)\\ H r + 2 

for any 2r > n. 

To complete the estimation, one observes that 

/ |?2(0,y)|dy< / \(j)(x,y)\dxdy = 

\y\>M \y\>M 

= I \fc*,y)\ (\i+V(\i+y\r (I 1 + x \y^ + y\) rdxd v < 

\y\>M 

(171) 

< C U\\l\I I (ll+xl yJ ill+yl ^ dxdy < 
y M>a/ 

So by setting M = , s\ = n + 1 in equation (|168|) . and applying the 
estimations ()169|1 . (|170j) . the result follows. 
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